Laser cooling of a diatomic molecule 
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It has been roughly three decades since laser cooling techniques produced ultracold atoms [J- 
[^, leading to rapid advances in a vast array of fields. Unfortunately laser cooling has not yet 
been extended to molecules because of their complex internal structure. However, this complexity 
makes molecules potentially useful for a wide range of applications For example, heteronuclear 
molecules possess permanent electric dipole moments which lead to long-range, tunable, anisotropic 
dipole-dipole interactions. The combination of the dipole-dipole interaction and the precise control 
over molecular degrees of freedom possible at ultracold temperatures make ultracold molecules at- 
tractive candidates for use in quantum simulation of condensed matter systems [3] and quantum 
computation la| . Also ultracold molecules could provide unique opportunities for studying chemical 
dynamics and for tests of fundamental symmetries [gl-ITll]. Here we experimentally demonstrate 
laser cooling of the polar molecule strontium monofluoride (SrF). Using an optical cycling scheme 
requiring only three lasers T^, we have observed both Sisyphus and Doppler cooling forces which 
have substantially reduced the transverse temperature of a SrF molecular beam. Currently the only 
technique for producing ultracold molecules is by binding together ultracold alkali atoms through 
Feshbach resonance or photoassociation [ij. By contrast, different proposed applications for 
ultracold molecules require a variety of molecular energy-level structures (e.g. unpaired electronic 
spin [sl. [gl. [Til. fisll ■ Omega doublets etc). Our method provides a new route to ultracold temper- 
atures for molecules. In particular it bridges the gap between ultracold temperatures and the ~ 1 
K temperatures attainable with directly cooled molecules (e.g. with cryogenic buffer gas cooling 
[TtI or decelerated supersonic beams [3). Ultimately our technique should enable the production 
of large samples of molecules at ultracold temperatures for species that are chemically distinct from 
bialkalis. 



Laser cooling of atoms has enabled unprecedented ac- 
cess to ultracold temperatures. The power of laser cool- 
ing generally derives from the ability of certain atoms 
to continuously scatter photons from a laser. For in- 
stance, Doppler laser cooling relies on small but repet- 
itive momentum kicks resulting from the absorption of 
red-detuned photons counterpropagating to the motion 
of an atom. Doppler cooling of an atom of mass m > 20 
amu with visible light from room temperature to ultra- 
cold temperatures requires > 10^ photon scatters. To 
scatter this many photons, an atom must have a closed 
cycling transition in which each photon absorption is al- 
ways followed by spontaneous decay back to the initial 
state. 

Unfortunately there are no completely closed transi- 
tions in any real physical systems, and inevitably spon- 
taneous decays to other states occur, usually before ul- 
tracold temperatures can be reached. Each additional 
populated level requires a "repump" laser to return pop- 
ulation back into the main cycle, so that photon scat- 
tering can continue. Cycling transitions requiring one 
or two "repump" lasers are common in atomic systems, 
but they are quite difficult to find in molecules because 
of their vibrational and rotational degrees of freedom. 
Control over vibrational states is particularly problem- 
atic because there is no strict selection rule governing the 
branching ratios for decay of an excited electronic state 
into different vibrational levels. Instead these branch- 
ing ratios are governed by the molecule's Franck-Condon 
factors (FCFs), which describe the overlap of the vibra- 



tional wavefunctions for different electronic states. For a 
typical molecule, the probability to return to the original 
vibrational level after 10^ photon scatters is extremely 
small. Furthermore, decay from a single excited state 
can populate up to three rotational levels per vibrational 
level, since rotational selection rules generally only re- 
quire that AA^ = 0,±1. (We use N, v, F, and M 
as the rotational, vibrational, total angular momentum, 
and Zeeman quantum numbers respectively.) Because 
each substantially populated level requires an individu- 
ally tunable "repump" laser, laser cooling of a molecule 
can easily require so many lasers as to be experimentally 
challenging. 

Here we experimentally demonstrate direct laser cool- 
ing of a diatomic molecule by reducing the transverse ve- 
locity spread of a cryogenic beam of strontium monoflu- 
oride (SrF) in ID. Cryogenic buffer gas beam sources 
produce highly directional beams with large fluxes for 
a variety of molecules, as described elsewhere [Hi [20| . 
Our scheme for eliminating the rotational and vibrational 
branching, necessary to ensure optical cycling, has been 
described previously [l^ . We briefly recount the main 
points in the context of the experiment reported here. 
We use the X^T,'^ ^ A'^Ili/2 electronic transition of SrF 
for cycling. Use of the first excited state A'^Ili/2 ensures 
that no other electronic states can be populated by spon- 
taneous decay [2lj. The state Jias a large spon- 
taneous decay rate, F = 27r x 7 MHz which enables 
the application of strong optical forces. We have chosen 
SrF primarily because its favorable FCF's dictate that 
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only three vibrational levels will be significantly popu- 
lated after 10^ photon scatters 23|, in principle more 
than sufficient for stopping molecules in our cryogenic 
beam. Rotational branching is eliminated by driving an 
= 1 — > A^' = type transition 24] , where the prime in- 
dicates the excited state. Transitions of this type lead to 
optical pumping into dark ground-state Zeeman sublevels 
not excited by the laser [25] . For example, in this system 
the F = 2 ]M] = 2 sublevels are dark when driven by 
linear laser polarization. We eliminate these dark states 
by applying a magnetic field oriented at an angle b with 
respect to the fixed linear laser polarization, forcing the 
dark states to Larmor precess into bright states. Finally, 
radiofrequency sidebands on the lasers address all ground 
state hyperfine (HFS) and spin-rotation (SR) substruc- 
ture. With this scheme only three cooling lasers (one 
main Ago = 663.3 nm pump laser and two Aio = 686.0 
nm and A21 — 685.4 nm vibrational repump lasers) are 
required as shown in Fig. 1. (Here Xij is the wavelength 
of the X{v = i) A{v' = j) transitions.) 

Use of the N = 1 ^ N' ^ type transition leads 
to a lower photon scattering rate than a traditional two- 
level atomic cycling transition 12|, |24| . To overcome this 
problem we use an elongated transverse cooling region 
(see Supplementary Information). The SrF beam is in- 
tersected by the three cooling laser beams at nearly right 
angles. The laser beams are reflected back and forth at 
a slight angle so that they intersect the SrF beam ~ 75 
times in the 15 cm long cooling region. At the end of 
the cooling region the beams are nearly retroreflected, 
resulting in the formation of standing waves. The com- 
bination of magnetic field remixing of Zeeman sublevels 
and standing waves can lead lead to Sisyphus forces in 
addition to Dop pler forces as has been observed in atomic 
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10 cm downstream from the end of the 
cooling region, laser induced fluorescence (LIF) is imaged 
to obtain the spatial distribution of the molecular beam. 
This spatial distribution maps onto the velocity distri- 
bution of the molecules with a resolution of ^ 1 m/s, 
so from such images we can extract information about 
the velocity-dependent forces applied to the molecules, 
as well as the beam's transverse temperature T. 

We find there are two cooling regimes with qualita- 
tively different features that depend critically on the 
magnitude of the applied B field. In Figs. 2a and 2b 
we show data with B = 5 G and 0.6 G respectively that 
are characteristic of these regimes. In both regimes we 
find that 9b is unimportant HOb 7^ 0, 90°. In all cases we 
observe that the total integrated LIF signal is constant to 
within the experimental reproducibility ('^ 5%), and so 
changes in the molecular spatial distribution accurately 
reflect changes in the velocity distribution of the molecu- 
lar beam. The top and bottom panels of Fig. 2a and 2b 
are representative molecular beam images for two differ- 
ent main pump laser carrier frequencies, Sp = ±1.5r. In 
this system there is not a single well-defined value of the 
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FIG. 1: Energy level structure in SrF. a) Relevant electronic 
and vibrational structure. Solid upward lines indicate laser- 
driven transitions in the experiment at wavelengths A„ „/. 
Wavy lines indicate spontaneous decays from the A{v = 0) 
state (solid) and the A{v = 1) state (dashed) with calcu- 
lated FCFs b) Relevant rotational energy levels, split- 
tings [26[ , and transitions (vertical arrows) in the SrF cycling 
scheme. SrF has an unpaired electron spin 5 = 1/2 which 
sphts the X(iV = 1) level into J = N ±S levels through SR 
interaction, and ^^F has nuclear spin 1 — 1/2 which splits 
J into F = J ± 1/2 through HFS interactions. HFS is not 
resolved in the A state |27|. c) Laser (vertical lines) and 
molecular (curved) spectra for addressing all HFS ground 
state sublevels. Molecular lines are shown with a power- 
broadened linewidth of l.SF and experimentally determined 
line strength. Electro-optic modulation of the laser frequency 
with index m = 2.6 and frequency /mod ~ 43 MHz produces 
four sidebands of nearly equal strength that nearly match the 
HFS splitting. The laser spectrum shown has fmod = 46.4 



MHz and Sr, 



-1.5r, which are experimentally determined 



to be optimal Aqo laser parameters for Doppler cooling. The 
laser carrier frequency corresponding to Sp = is defined 
experimentally by the frequency (thick central line in laser 
spectrum) which produces maximal LIF, and 5p corresponds 
to the detuning of the carrier laser frequency from this value. 



detuning for the pump and repump lasers because each 
laser contains several sideband frequencies that each in- 
teract with multiple transitions between ground and ex- 
cited states. We define Sp(^j.) = experimentally by deter- 
mining the pump (repump) laser carrier frequency which 
produces maximal LIF. For 6p = — (+)1.5F the detun- 
ing of the nearest sideband from each transition ranges 
from to ^ — (-f )1.5F, indicating that the sign of Sp cor- 
responds to net average red (blue) detuning. (Here and 
throughout the rest of the paper both vibrational repump 
lasers have modulation frequencies /^^^ = 43 MHz and 
6r = 0. The main pump laser has modulation frequency 
fmod ^ listed in the captions.) 
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For a red detuned main pump laser, 6p = — l.SF, and 
B = 5 G we observe significant narrowing of the molecu- 
lar beam and enhancement of molecules with low trans- 
verse velocity, Vt , as shown in the bottom panel of Fig. 
2a. This corresponds to a reduction in the spread of vt , 
and is a clear signature of Doppler cooling. Also it is ev- 
ident that the entire molecular beam experiences cooling 
forces, indicating that the cooling force is significant for 
all Vt in the molecular beam. The molecular beam is con- 
strained by collimating apertures to have Vt < i'mTx" — 4 
m/s. For a blue detuned main pump laser, Sp = -I-1.5F, 
and -B = 5 G we observe depletion of low molecules and 
broadening of the molecular beam as shown in the top 
panel of Fig. 2a, as expected for Doppler heating. Un- 
der these conditions, there is also a small but noticeable 
sharp feature in the center of the molecular beam, indi- 
cating some residual cooling of the remaining molecules 
with low Vt- 

In Fig. 2b we show data characteristic of a small ap- 
plied B field of 0.6 G, which are strikingly different from 
the data in Fig. 2a. The most significant difference be- 
tween Figs. 2a and 2b is that cooling occurs for detun- 
ings of opposite sign. For a red detuned main pump 
laser, Sp = — l.SF, we observe (Fig. 2b, lower panel) 
two relatively sharp peaks, neither of which is centered 
around Vt = 0- This indicates the heating of molecules 
with low \vt \ and accumulation around two stable veloc- 
ity points with vt ^ 0. For a blue detuned main pump 
laser, 6p > —l.SF, we observe (Fig. 2b, lower panel) a 
sharp central spike and a large enhancement of low vt 
molecules. This feature results from the strong cooling 
of molecules over a small range of Vt around zero. Mean- 
while molecules with larger \vt \ experience small heating 
forces, resulting in a very slight enhancement at large 
\vt\. 

A complete characterization of the detailed cooling 
forces responsible for these observations would require 
the solution of the optical Bloch equations for this sys- 
tem. All relevant quantities (detunings, Larmor frequen- 
cies, Rabi frequencies, F) are the same within factors 
of order unity, so the full 44 level system, driven by 
twelve laser frequencies, each interacting with multiple 
levels, must be solved. Such a calculation is beyond 
the scope of this work; however, the main features we 
have observed are common to any system with magnet- 
ically remixed dark sublevels driven by a standing wave 
2^33|- The simplest such system is an _F = 1 
transition driven by a single laser frequency, and it pro- 
vides substantial insight into our observations. A lin- 
early polarized, blue (red) detuned laser of wavelength A 
produces an AC Stark shift which attracts (repels) the 
i^' = 0,A/ = and F' = 1,M' = levels, while the 
F — 1, M — ±1 levels are unperturbed by the light field 
to first order, as shown schematically in the Supplemen- 
tal Information. In a standing wave, the F = 1,M = 
level undergoes a spatially periodic AC Stark shift with 



period A/2. For blue (red) detuning, molecules ride up 
a potential hill in this level, losing (gaining) kinetic en- 
ergy, before they are pumped at a rate 7p < F into the 
F — 1,M = ±1 sublevels which are dark. If a i3-field 
is applied a,t 9b ^ 0°,90° then the molecules precess 
from the dark A/ = ±1 sublevels back into M = at 
the nodes of the standing wave at a rate ujb- The Sisy- 



phus force is maximized when vt 



f7p, and 



ojSjs ^ (-ggg Supplemental Information), vf}^^ defines 
the effective velocity range of the Sisyphus force. Larger 
values of produce a much smaller Sisyphus force be- 
cause the molecules are pumped back and forth between 
bright and dark states at random points in the stand- 
ing wave. Because the molecules spend more time in the 
bright states the photon scattering rate is higher, and 
the Doppler force is larger. Simple arguments suggest 
that ujg°^ ~ F for maximum Doppler cooling forces, and 
that Doppler forces occur over a larger range of velocities 
Vt — v^°^^ = ^=4^ ni/s (see Supplemental Information). 

This qualitative discussion provides substantial insight 
into our observations. At low B fields we expect to ob- 
serve Sisyphus forces, which are characterized by cooling 
(heating) for blue (red) detuning for molecules with low 
Vt as observed in Fig. 2b. At higher B fields we ex- 
pect to observe Doppler forces, characterized by cooling 
(heating) for red (blue) detuning. Since vf/^™ ~ "^a/Ix^ 
Doppler cooling forces should affect the entire molecular 
beam as observed in the bottom panel of Fig. 2a. In 
the Supplemental Information we provide the argument 
for the estimate 7p ^ ^, which yields vf}^^ ^ 0.2 m/s. 
This is much smaller than wf/J'J" — ~ 4 m/s, and 

is consistent with our observations. 

Of course neither regime can be characterized purely 
by either Sisyphus or Doppler forces. In the moderate 
B regime residual Sisyphus forces lead to a slight addi- 
tional broadening (narrowing) in the low-velocity part of 
the distribution for red (blue) detuning. This residual 
narrowing is clear in the top panel of Fig. 2a, while the 
broadening is too small to observe in the bottom panel 
of Fig. 2a. In the low B regime residual Doppler forces 
result in small heating (cooling) over a large range of 
velocities for blue (red) detuning. This gives rise to non- 
zero unstable (stable) velocities vt — iuo, where the net 
force is zero and population is depleted (accumulates). 
These stable points vt = ±vo are clear in the bottom 
panel of Fig. 2b. In the top panel of Fig. 2b the un- 
stable velocities are clearly depleted, and there is some 
residual Doppler broadening in the wings of the molecu- 
lar beam. Although the arguments presented here derive 
from an F = 1 — > _F' = example system, they are com- 
mon to any F —i' F' < F system. Furthermore the effects 
described here have been observed in several such atomic 
systems [^ISl- 

To clearly illustrate the B field dependence of the cool- 
ing force, we show the magnetic field dependence of the 



4 




cooling forces for a red detuned pump laser, Sp = — l.SF, 
in Fig. 3a. As shown for very small B fields, the width 
of the molecular beam increases due to Sisyphus heat- 
ing effects. This increase is followed by a sharp decrease 
in the molecular beam width corresponding to Doppler 
cooling. The width of the molecular beam is quite insen- 
sitive to the magnetic field over a range of intermediate 
magnetic field amplitudes between 2 and 6 G. Finally at 
magnetic fields higher than 6 G the Doppler forces are 
reduced because the magnetic field artificially broadens 
the transitions, resulting in lower scattering rates, and 
lower Doppler forces. We estimate values for in the 
Supplemental Information which yield Boop ~ 5 G and 
Bsis 0.2 G for the maximum Sisyphus and Doppler 
forces, in reasonable agreement with our observations. 
Once again, detailed comparison would require the full 
solution of the optical Bloch equations for our system. 

As an example of the complex features present in this 
system, we show the frequency dependence of the width 
of the SrF beam under Doppler force-dominated condi- 
tions of B = 5 G, for various pump laser detunings Sp in 
Fig. 3b. The frequency dependence of the width is sub- 
stantially more complicated than that of a typical 2-level 
system. As shown the force oscillates many times be- 
tween heating and cooling for —250 MHz < 6p < 250 
MHz. However, this complicated structure is amenable 
to a simple interpretation. As the pump laser frequency 
is varied, the nearest laser frequency to each molecular 
transition oscillates between red-detuned (cooling) and 
blue-detuned (heating) with a variation that is nearly 
periodic in the sideband frequency. The frequency de- 
pendence in the Sisyphus regime, shown in Fig. 3c, is 
somewhat more complicated, but for small detunings —50 
MHz< 5p < 50 MHz the Sisyphus force has the opposite 
sign as the Doppler forces as expected. In the Supple- 
mental Information we also show the power dependence 
of the Sisyphus and Doppler cooling forces which are both 
in reasonable agreement with expectations. 

Finally, we discuss our determination of the tempera- 
ture T of the molecules after they are cooled. The unper- 
turbed beam is constrained by collimating apertures to 
have To = 50 mK. The long interaction region prevents a 
precise determination of T for the cooled beam, because 
the beam continues to expand as it experiences an im- 
perfectly known distribution of cooling forces throughout 
the interaction region. To estimate T, we calculate the 
molecular beam profile using a Monte Carlo simulation of 
classical particles subjected to the qualitatively expected 
force vs. velocity profile (see Supplemental Information). 
Using these simulations we find T ~ Tsis — 300/iK for 
the Sisyphus regime in the top panel of Fig. 2b. We also 
estimate a conservative upper limit on the temperature 
of T < T|ff ^ = 5 mK. 

For the Doppler regime, we find that T ~ Tbop = 5 
mK, and T < T]!}^^ = 15 mK. These values of T for 
the Doppler regime are consistent with the final temper- 



FIG. 2: Laser cooling of SrF. LIF in the probe region with- 
out cooling lasers in the interaction region (black curves), 
with cooling lasers and main pump laser red-detuned by 
Sp — — 1.5r (red), and with cooling lasers and main pump 
laser blue-detuned by Sp = -hl.SF (blue) for (a) = 46.4, 

B = 5 G, 0s = 60°, and (b) ff^^^ = 43.2 MHz, B = 0.6 
G, 9b ~ 30°. In (a) cooling (heating) of the beam is ob- 
served for red (blue) detuning; both are in accordance with 
expectations for Doppler forces. In (b) cooling (heating) is 
observed for blue (red) detuning, as is expected for Sisyphus 
forces. In all cases the total integrated signal is the same to 
within 5%. The blue-detuned curves systematically have the 
lowest total integrated signal, likely because some molecules 
have been been pushed outside the imaging region by the 
Doppler heating force. We observe qualitatively similar be- 
havior in both cases for pump laser modulation frequencies 
42 < /^o^ < 47 MHz. However the /^^^j values shown here 
were chosen because they produced the clearest Doppler and 
Sisyphus effects. This plot shows raw data, with no rescaling 
applied. 



ature expected if the molecules are subjected to N^c ~ 
500 — 1000 photon scatters. This value of Ngc agrees with 
expectations based on the previously observed scattering 
rate for this system and roughly known interaction time. 
Importantly, the total integrated signal of the Doppler- 
cooled beam and the unperturbed beam are the same to 
within the level of experimental reproducibility (« 5%), 
indicating that our cycling scheme is highly closed. 

Our results have immediate implications for a num- 
ber of future experiments. For example, the ID cool- 
ing and optical cycling demonstrated here could dra- 
matically improve the statistical sensitivity of searches 
for electron electric dipole and nuclear anapole moments 
llH, by providing more collimated molecular beams 
and enhanced detection efficiency in these experiments. 
In addition, the combination of ID cooling and a highly 
closed cycling transition opens the door to laser cooling 
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of a wide variety of new physical phenomena. 



FIG. 3: Magnetic field and frequency dependence of the cool- 
ing forces, a) Magnetic field dependence of the cooling force. 
The width of the molecular beam is plotted for f^^^ ~ 46.4 
MHz, Sp = -1.5r, and various B fields at Ob = 30°. For 
B < 2 G the beam is broader due to Sisyphus forces. At ap- 
proximately B = 2 G we see a sharp transition between Sisy- 
phus dominated heating forces, and Doppler dominated cool- 
ing forces. For B > 6 G the width increases as described in the 
main text. Frequency dependence of the cooling forces under 
b) Doppler {B — 5 G) and c) Sisyphus {B — 0.6 G) dominated 
conditions. The solid lines between data points are drawn as 
a guide to the eye. The vertical dashed lines are spaced by 
the sideband frequency of a) 46.4 MHz and b) 43.2 MHz, and 
illustrate the dependence of the oscillations on sideband fre- 
quency. Describing the molecular beam simply by a width 
is potentially problematic in the Sisyphus dominated regime 
where more than one velocity class are present. Nonetheless 
the total width of the beam for the red (blue) detuned pump 
laser is larger (smaller) than the unperturbed beam width, 
and thus the width provides a reasonable description of the 
molecular beam. The asymmetry around 5p = in b) and c) 
underscores the complex nature of this system and is likely 
representative of the fact that every sideband does not have 
the same detuning from the nearest transition. Therefore, 
the detuning of each sideband from the nearest transition is 
not zero for Sp = 0. Furthermore, the absolute values of the 
individual sideband detunings are not the same for +Sp and 
—5p. 



of molecules in 3D. Given the calculated FCF's, a large 
fraction of molecules should scatter the ~ 40, 000 photons 
necessary to bring a beam of SrF to a stop, and subse- 
quently load the molecules into a trap. Furthermore the 
experimentally determined loss rate of molecules in this 
system is < 5% for Nsc ~1000, implying that > 10% 
can be brought to rest, given sufficient interaction time. 
The laser cooling techniques presented here are limited, 
from a practical standpoint, to those molecules which 
have closed electronic transitions with diagonal FCFs and 
therefore require relatively few lasers. For this reason 
these techniques are applicable only to a small fraction 
of diatomic molecules. However because the set of di- 
atomic molecules is very large, this subset contains a sig- 
nificant number of molecules [23|. We are aware of per- 
haps a dozen diatomic molecules with a wide range of 
internal structures that appear amenable to laser cooling 
with similar methods. Laser cooling such molecules to ul- 
tracold temperatures would open the door to the study 
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SUPPLEMENTARY INFORMATION 
1 Experimental setup 

A schematic of our experiment is shown in Supplemen- 
tary Fig. |4l The source of SrF in our experiment is a 
(2.5x2.5x2.5 cm'^) cryogenic buffer gas cell. The cell is 
attached to the 4K surface of a liquid helium-filled dewar 
inside a vacuum chamber held at 5 x 10^^ Torr. 4K 
helium buffer gas flows continuously into the cell through 
an inlet in one face of the cell. We typically flow 5 seem 
(standard cubic centimeters per minute) of helium gas 
into the cell which leads to a steady state density of 
~ 10^^ cm^'^. We introduce SrF into the cell by ablating a 
solid SrF2 target formed by pressing granular SrF2 pow- 
der into a disk. The target is ablated using 5 ns, ^20 mJ 
pulses of 1064 nm light from a Q-switched Nd:YAG laser 
typically operating at 1 Hz repetition rate. A mixture of 
particles including neutral SrF is ejected from the target 
and thermalizes to 4K after sufficient collisions with the 
helium buffer gas. A 3 mm diameter exit aperture on the 
cell wall directly opposite the helium inlet is used to form 
the molecular beam. The beam consists of SrF, helium 
and the other particles created by ablation. The SrF in 
the beam has a measured forward velocity of U| | = 200±30 
m/s with a spread of Au|| «40 m/s. The beam has ~ 10^ 
molecules in the X(i; = 0, 7V = 1) state per shot. 

2.5 cm downstream from the cell, a coconut-charcoal 
covered copper plate with a 6.25 mm diameter hole colli- 
mates the SrF beam while also limiting the flow of helium 
into the rest of the apparatus. The beam then exits the 
cryostat where it is further coUimated by a 1.5 x 1.5 mm^ 
aperture, resulting in an RMS transverse velocity spread 
yRMS — i Q ni/s. 3 cm after the coUimating aperture, 
the beam enters the cooling region, which is defined by 
15 cm long windows. The cooling laser beams intersect 
the SrF beam at nearly right angles. The linearly polar- 




2.5 11 14 29 39 cm 

FIG. 4: Schematic of the experimental apparatus. 



ized laser beams are reflected back and forth using 20 cm 
long mirrors located outside the vacuum chamber. The 
windows are mounted at Brewster's angle to allow the 
laser beams to pass multiple times through the cooling 
region with minimal attenuation. We apply -B-fields us- 
ing two pairs of rectangular coils mounted outside the 
cooling region. LIF can be monitored at any position in 
the cooling region using a photomultiplier tube. 

10 cm after the end of the cooling region, the SrF beam 
enters the probe region where LIF from the X{v — 0,N — 
1) -H> A{v' = 0,J' = 1/2) transition is imaged onto an 
intensified CCD camera. A "cleanup" Aio laser beam 
intersects the SrF beam between the cooling and probe 
regions to return any residual population in the X(w — 
1) state back to the X{v — 0) state. Both probe and 
"cleanup" beams are retroreflected to eliminate artificial 
Doppler shifts. 

2 Capture velocity and optimum magnetic Held for 
Sisyphus and Doppler forces 

In this section we provide estimations of the capture 
velocities and optimum magnetic fields for Sisyphus and 
Doppler forces, which should be accurate within a factor 
of ^ 2. Once optical cycling is achieved in a system, 
the application of laser beams in opposite directions can 
lead to substantial cooling forces. The simplest of these 
is the Doppler force, where the Doppler shift of a moving 
molecule brings it closer to or further from resonance 
with the laser. Molecules closer to resonance scatter more 
photons, leading to a velocity-dependent force, Fuop oc 
vt- The Doppler force affects molecules with velocities 
such that the Doppler shift, ujd = kvt, is not greater 
than the natural linewidth F, where k = 27r/Aoo- This 
leads to an effective velocity range for the Doppler force 
given by: \vt\ < v^°ax — ^/^ — ^ ^/s- Molecules with 
larger velocities than w^"*^ experience a reduced force. 

In systems with dark Zeeman sublevels, such as the 
SrF and F = 1 F' = systems presented in the main 
text, optical cycling ceases as soon as the molecules are 
pumped into these dark states. If a magnetic field is ap- 
plied then the dark states can Larmor precess into bright 
states and optical cycling can continue. The applied mag- 
netic field must be large enough that the molecules evolve 
out of the dark states quickly, but it must not be so large 
that they precess out of the bright states before they 
can absorb a photon. An estimate for the optimal B 
field for Doppler cooling, Bjjop, can be made by equat- 
ing the Larmor precession frequency, lob ~ BuBgp, and 
the linewidth of the transition F, where /is is the Bohr 
magneton (the magnetic moment of SrF) and (7f is the 
Lande g factor for each hyperfine level. In SrF, gp 1 
for the hyperfine levels. We then have B]jop ^ — = 5 
G. ^'^ 

As described in the main text, these types of systems 
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also give rise to Sisyphus forces if they are subjected to a 
standing wave. The basic mechanism responsible for the 
Sisyphus effects is shown schematically in Supplemen- 
tary Fig. [S] Sisyphus forces are maximal for molecules 
that travel a distance of A/4 in the time it takes for the 
molecules to be pumped into the dark states, — . In 
this case the laser field extracts the maximum kinetic 
energy from the molecules as they ride up the entire po- 
tential hill. This leads to an effective velocity range of 
= vf^lx ~ i7p fo"^ Sisyphus forces. To estimate vf}^^ 
we need to determine 7p. In the SrF system presented 
in the main text, only the F = 2, \M\ = 2 sublevels 
are dark in general, and a simple counting argument can 
give an estimate for 7^ in this case. If all transitions 
are completely saturated, then the molecules spend equal 
time in the 24 v — 1 and v = ground states, and the 
4 excited states. The maximum photon scattering rate 
is then 711^^: = ^ ^- According to the calculated 

FCF's shown in the main text, 98% of the spontaneously 
emitted photons from the excited state result in decays 
to the V = ground state sublevels. To first order, 1/6 
of spontaneous emission events result in the population 
of the \M\ = 2 sublevels. We then estimate 7p ~ 
or 7p ^ and we find that vf}^^ 0.2m/s. 



a) 

F'=0 



The previous discussion gives the optimal vt for Sisy- 
phus forces; however, the maximum force can only be 
achieved for over a small range of B fields. If B is too 
small or too large, then the molecules precess back and 
forth between bright and dark states at random points 
in the standing wave. The net result is that after sev- 
eral Sisyphus cycles, the molecules have neither gained 
nor lost a significant amount of energy. If on the other 



hand ujb — 7p then molecules with vt 



Sis 



travel A/4 



in the time it takes to precess out of the dark states. 
At this point the molecules can repeat the Sisyphus pro- 
cess shown in Supplementary Fig. [SJj. Since only the 
\M\ =2 sublevels are dark, we only need to consider 
their precession rate to estimate the optimal B field for 
Sisyphus cooling. For F = 2, gp — 1/2, and we obtain 



F = l 




m=-l m=0 m=l 



-15 -10 -5 5 10 15 
Velocity(m/s) 

FIG. 5: Schematic illustrating the origin of Sisyphus force in 
an F = 1 F' = Q system, (a) Zero field energy levels 
(solid) and AC Stark shifts (dashed) in the presence of a blue 
detuned, linearly polarized laser field, (b) Energy levels and 
motion of molecules in a blue detuned standing wave. As ex- 
plained in the text, the molecules effectively ride continuously 
up potential hills, losing kinetic energy. For red detuning the 
picture is reversed, with molecules gaining kinetic energy as 
they ride down potential hills. For optimum Sisyphus forces, 
the molecules must traverse from node to antinode in the op- 
tical pumping time, l/7p, and the molecules must traverse 
from antinode to node in the magnetic field remixing time, 
1 / u)B ■ (c) Results of a semi-classical calculation of the average 
force on a molecule in a standing wave with Rabi frequency 
Q,R — r, and red detuning S = —T for low (bottom), interme- 
diate (middle), and high (top) field remixing rates ujb- The 
curves are offset by 0.4 hk/F for clarity. For small ujb we ob- 
serve Sisyphus heating at small velocities and small Doppler 
cooling at large velocities, with two stable points at vt — ±2.5 
m/s. For higher ujb observe almost pure Doppler cooling 
forces for all velocities. For the highest ujb the Doppler force 
decreases because the B-field broadens the transitions result- 
ing in lower photon scattering rates. For a blue detuning of 
6 = +r, the force is the same magnitude, but reversed in sign. 
The calculation then predicts Sisyphus cooling for small ujb 
and Doppler heating for higher ujb- 



3 Estimation of Temperature 

As mentioned in the main text, we use a Monte Carlo 
simulation of classical particles subjected to the quali- 
tatively expected Sisyphus and Doppler force versus ve- 
locity curves. To obtain qualitative estimates of these 
force curves we have solved the optical Bloch equations 
for an F = 1 — > F' = system in a standing wave and 
magnetic field. We then compute the average force over 
one wavelength. Typical results of such calculations are 
shown in Supplementary Fig. [5t, and show good agree- 
ment with the qualitative discussion in the main text. 
Specifically, in the calculation we see that there are three 



distinct magnetic field regimes. At low B fields the cal- 
culation shows strong Sisyphus forces. For low B fields, 
blue (red) detuning produces cooling (heating) forces for 
\vt\ < wffaa; — 0.4 m/s. We also see that for \vt\ > 2 
m/s there are residual Doppler forces that are of oppo- 
site sign as the Sisyphus force. This gives rise to non-zero 
unstable (stable) velocities Vt — ±2 m/s, where the net 
force is zero and population is diminished (accumulates). 
We also see that Sisyphus forces are not substantial for 
\vt \ > wf/aa;, SO vf}l^ indicates the velocity extent of the 
Sisyphus force. 

At intermediate B fields the Sisyphus forces are sub- 
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FIG. 6: Results of a Monte Carlo simulation of the SrF beam 
subjected to Sisyphus type forces, a) Experimental (noisy) 
and Monte Carlo simulation (clean) molecular signal for the 
unperturbed (black), Sisyphus cooled (blue), and Sisyphus 
heated (red) molecular beam. The curves are offset by one 
unit for clarity. The total area of the simulated molecular 
beam curves are normalized to their respective experimen- 
tal curve. The difference in area between the experimental 
Sisyphus cooled beam and the unperturbed beam is ~ 5%. 
b) Velocity distribution (solid) and Gaussian fit (dashed) for 
the unperturbed (black) and Sisyphus-cooled (blue) molecu- 
lar beam resulting from the Monte Carlo simulation of the 
molecular beam. From the Gaussian fit we extract the trans- 
verse temperature T of the molecular beam. 



stantially reduced, and instead the force curve is indica- 
tive of Doppler forces. In contrast to the Sisyphus mech- 
anism, red detuning produces Doppler coohng, while 
blue detuning produces Doppler heating. Furthermore, 
the Doppler cooling/heating forces extend over much 
larger transverse velocities than the Sisyphus forces. The 
Doppler forces have a broad maximum value around 
t)4 ~ 4 m/s, indicating that ^ 4 m/s as expected. 

At this intermediate B field, there still remains a small 
residual Sisyphus force for small Vt that has opposite sign 
to the Doppler force. At the highest B fields the calcu- 
lation shows reduced Doppler forces because the mag- 
netic field artificially broadens the transitions, resulting 
in lower scattering rates, and lower Doppler forces. 

We have used the general shape of the force curve de- 
rived from the calculation in our simulations of the molec- 
ular beam after exposure to Sisyphus and Doppler type 
cooling conditions. In these simulations we have used 
the value of v^^^ obtained from the F ^ 1 ^ F' = 
calculation. In the calculation the fraction of dark states 
(in the absence of the remixing S-field) is larger than in 
the real SrF system; hence we take the value of v^^^ ob- 
served in the calculations as an upper limit on the value 
in our experimental system. To estimate the temperature 
of the Sisyphus cooled beam, we assume that the Sisy- 
phus cooling occurs only over the last third of the cooling 



region (where the standing waves are most pronounced). 
We then adjust the magnitude of the force by an overall 
factor l3 until the simulation matches the experimental 
LIF profile. 

Examples of these simulations are shown in Supple- 
mentary Fig. under Sisyphus conditions, and show 
good agreement with the experimental results. The final 
temperature estimated in the case of Sisyphus cooling is 
primarily dependent on v^^^. v^^^ is proportional to 7p, 
and the temperature extracted is higher for larger values 
of this velocity. To obtain a conservative estimate of T 
we have used the value of v^^^ from the F — 1 —i' F' — 
calculation. This value is likely larger than in the actual 
SrF system as shown by our estimation, so we expect 
that this calculation tends to overestimate T. Using this 
simulation we find T ~ Tsis — 300/iK. However, because 
the full SrF system was not used to calculate the force 
used in the simulation, we derive a bound for the max- 
imum value of T by making two overly conservative as- 
sumptions. First, we assume that the cooling force only 
occurs over a 1 cm length at the beginning of the interac- 
tion region. Second, we allow v^^^ to be larger than the 



value found from either the F 



^ cap 
1 



F' = calculation 



or our estimation. We then find the largest value of v^l^ 
which still replicates the experimental data under varia- 
tion of /3. Under these very conservative assumptions we 

This value is 



estimate T|f/^ 



= 5 mK for v^^^^ = 0.8 m/s. 

Sis 



significantly larger than either v^^^ = 0.2 m/s from our 
estimation, or v^l^ = 0.4 m/s from the F = 1 ^ F' = Q 
calculation. 

In the Doppler regime, we assume a uniform cooling 
force F = —av over the whole interaction region and 
molecular beam velocity distribution. The overall mag- 
nitude of the force is adjusted until the simulated molec- 
ular beam width matches the width of the experimental 
LIF profile. With this method wc find T ~ Toop = 5 
mK. We also estimate a conservative upper limit by as- 
suming this cooling force is only applied over a 1 cm 
length in the beginning of the cooling region; from this 
we find T < T^^^ — 15 mK. If we use the more realistic 
force vs. velocity curve shown in Supplementary Fig. [S]c 
rather than F = —av we obtain similar temperature es- 
timates, and slightly better agreement between the simu- 
lation and the data. Specifically, inclusion of the residual 
Sisyphus force in the Doppler force curve reproduces the 
small sharp feature in the center of the molecular beam 
in Fig. 2a of the main text. 

We have also performed calculations for a system with 
lower states J = 3/2 and J — 1/2 excited to a single 
J' = 1/2 upper state driven by two laser frequencies. 
In this calculation we assume that each laser only drives 
population from the nearest J level. This calculation 
produced a very similar i?-field dependence and temper- 
ature estimate as the calculation based on the simpler 
system. 
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In Fig. [7] we have plotted the pump laser power depen- 
dence under both Sisyphus and Doppler cooling condi- 
tions with 6p = +1.5r (Sisyphus cooling) and 6p = — l.SF 
(Doppler cooling). As shown in Supplementary Fig. [7)d, 
the Doppler cooling force increases linearly with power 
until it becomes saturated. This type of behavior is ex- 
pected for Doppler cooling because at low powers the 
Doppler force varies as the photon scattering rate which 
is linear with power. At higher powers, the photon scat- 
tering rate saturates because it is limited by the spon- 
taneous emission rate, and larger excitation rates do not 
lead to larger photon scattering rates. 



FIG. 7: Molecular beam temperature under a) Sisyphus and 
b) Doppler conditions as the power in the main pump laser 
is varied. Doppler cooling saturates at a main laser power of 
approximately 30 mW, while the Sisyphus cooling continues 
to increase as the power is raised. At powers below 10 mW, 
Sisyphus cooling is no longer observed. 



4 Power dependence of Sisyphus and Doppler 
cooling 

In this section we present additional data illustrating 
the dependence of the cooling forces on the laser power. 



The Sisyphus cooling force, as shown in Fig. [7^, how- 
ever, does not appear to saturate. The energy extracted 
on each Sisyphus cycle is proportional to the AC Stark 
shift of the ground state. For Rabi frequencies smaller 
than the detuning, the AC Stark shift varies as the inten- 
sity of the laser, while for Rabi frequencies larger than 
the detuning the AC Stark shift varies as the electric 
field of the laser. The data shown here are not sufficient 
to draw firm conclusions regarding a linear or quadratic 
power dependence; however, it is apparent that Sisyphus 
cooling requires more power than Doppler cooling. 



